Hamiltonian description of self-consistent wave-particle dynamics in a periodic structure 
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The coupled dynamics of electrons and electromagnetic fields propagating in traveling wave tubes is 
expressed with a hamiltonian formulation. The field is represented with eigenfunctions adapted to Floquet 
boundary conditions along the tube axis, using the Gel'fand /3-transform. The electron hamiltonian is 
the standard one coupling the particles to the propagating fields. The dynamics conserves energy, and 
excludes self-acceleration. A complete hamiltonian formulation of the dynamics results from adding 
space charge effects by electrostatic action-at-a-distance coupling. 
PACS numbers: 84.40. Fe (microwave tubes) 
52.35.Fp (Plasma: electrostatic waves and oscillations) 
52.40.Mj (particle beam interaction in plasmas) 
52.20.Dq (particle orbits) 
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To model consistently the interaction between electrons 
and waves in devices such as the traveling wave tube, free 
electron laser or synchrotron, one is presently left with two 
options. The first one 0, B [HI] is to consider the flow of 
electrons as a distributed charge and current density cou- 
pled with the field through the Maxwell equations. Since 
it generates diverging singularities, the particle nature of 
electrons is intentionally overlooked until the question of 
determining the trajectories of the flow is raised. Here the 
only possibility is to switch back to a particle description 
in which the Lorentz force applies. This change of model 
for the flow precludes the description of the wave-electron 
system in hamiltonian form. One reason is that a procedure 
is needed to distribute the electron charge and current into 
a finite volume. This procedure, usually based on mes hing 
space, can only be arbitrary. The second option Jf| HI [HI 
is to consistently consider electrons as particles and to de- 
termine the field they radiate starting from the Lienard- 
Wiechert potentials. The problem is that it is not possible 
to calculate the reaction from the radiated field because all 
quantities are diverging at the electron position. Attempts 
to balance energy and momentum resorting to a posteri- 
ori forces such as the Abraham-Lorentz force are plagued 
with a number of difficulties [H, Ojl [l4 , 16 1- Among them 
the infinite rest mass of the electron, self- acceleration and 
acausality, are also incompatible with the existence of a 
well-posed hamiltonian. 

Hamiltonians are essential to consistently define energy 
and momentum. They also have great practical usefulness 
to find approximate solutions in complex systems or to con- 
trol errors in numerical integration schemes [0. So their 
absence in the case of the classical wave-electron inter- 
action is both theoretically and practically unsatisfactory. 



While the final solution to these problems may involve an 
(upgraded, regular) quantum electrodynamics, one could 
at least expect a consistent classical approximation, which 
would be the classical limit of its quantum counterpart. For 
instance, in the limit where the electron radiates a large 
number of photons, each of them having a small energy 
compared to its kinetic energy, which is the case for the 
aforementioned devices, one could expect that averaging 
the quantum theory over a large number of emitted pho- 
tons would provide a consistent classical limit. But such 
a theory, free of the above problems, remains to be found 
HO- To progress in this direction, we undertake the con- 
struction of a hamiltonian for the traveling wave tube. 

First of all, one can distinguish the space-charge field, 
which is the curl-free part of the fields, from the ra- 
diated field, which is the divergence-free part. In free 
space, the space-charge field between several electrons 
with positions derives from the hamiltonian 7-L sc = 
8^T e2 /ll r i — r j\\- Adopting this point of view elim- 
inates the problem of infinite rest mass, by eliminating al- 
together the electrostatic field degrees of freedom. So one 
is left with finding a hamiltonian for the radiated fields cou- 
pled with the electrons. This is the object of this work. 

The traveling wave tube is a waveguide, periodic along 
the x axis with period d. All equations are written in the 
waveguide reference frame. It is bounded laterally by per- 
fect metal (E x ei = 0, H ■ ei = 0) boundary condition. 
All periods, with shape Vn and volume j Vn |, form a domain 
Vz '■= U„ e z{r + nde x : r € V }- The waves propagat- 
ing in the tube are described with their electric and mag- 
netic fields, E(r, t), H(r, t), obeying the Maxwell equa- 
tions with sources inside the tube [Q]]. 

We decompose the radiated field on the function ba- 
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sis made of the propagating modes. These modes satisfy 
the lateral boundary condition and the Floquet condition, 
E(r + de x ) = e" i/3d E(r) for some (3 £ R. Inside V , 
they satisfy V • E S(3 = 0, V • H s/3 = 0, and 
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where the O s/ g are the real eigenvalues, fi Q and e the per- 
meability and permittivity of vacuum. For a given /3 £ 
[0, 2ir/d[, the set of eigen electric fields E S| g (s € N) 
is an orthogonal basis of the divergence-free subset of 
-ff (rot , Vo) (this is the Hilbert space of square-integrable 
fields on Vo which have a square-integrable curl). Any 
divergence-free field E /3 (r,t) on Vo can be decomposed 
on the E S(3 with coefficients V a p, 
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This decomposition can be continued to V% thanks to the 
Floquet condition, and the continuation field satisfies the 
Floquet condition as well (possibly with a Gibbs phe- 
nomenon, depending on the field regularity). 

An arbitrary function G in the propagating structure does 
not, in general, satisfy the Floquet condition. But it can be 
expanded on functions Gp satisfying the Floquet condition. 
Indeed, the Fourier series on the variable /3, with x as a 
parameter, 



G„(x,() 
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satisfies the Floquet condition with respect to x. The n-th 
Fourier coefficient in the series is given by 

»2tt 



i r* 

G(x + nd,t) = — / Gp(x,t) 

^ Jpd=0 



d((3d), (5) 



which for n = yields exactly the requested expansion. 
The transform ©, hereafter referred to as the /3-transform 
decomposes any function into a superposi- 
tion of functions satisfying the Floquet condition. It will 
be central in what follows. Relation d5]l is the inverse (3- 
transform. 

Using equation (J3j, we finally get 

E ( r >*) = ^E r V sP (t)E s0 (r)d((3d). (6) 



set- 



This decomposition holds for any free electric field in 
the structure, e.g. propagating modes, evanescent modes 
or superposition of these. The magnetic field expansion is 
written 
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with the factor i for later convenience. 

Our initial choice to consider independently the propa- 
gating field holds only if it is decoupled from the space- 
charge field. This latter adds the curl-free term — V(/> to 



the expression of E where </> satisfies the Poisson equa- 
tion A(f> = —p/e . The /3-transform commutes both with 
time and space derivatives, so Maxwell equations are valid 
for E^ and as well. Using the /3-transform of <j), the 
general form of the electric field in presence of charged 
particles is 



E0(r,t)=£V, /J (i)E, /) (r)-VV 



(8) 



Using Maxwell equations, the field decompositions (J8]) and 
(0, and the definition of eigenfields E S(3 and H S| a, we get 

— y~] IspQsp~Eisi3 = ^ V S/ 3E S(3 H — g^-, (9) 



at 



E 



(10) 



where J is the current density. We multiply the first ex- 
pression by the conjugate of the field E* //3 and the second 
expression by H*,o. Integrating the resulting expressions 
over the volume Vo yields vanishing terms for s ^ s' . 
Otherwise we exhibit the electric ~€q J Vq |E s ^ j 2 d 3 r and 

magnetic |/U J Vo |H s/ g| 2 d 3 r energies stored in one pe- 
riod for the corresponding propagation mode. Both en- 
ergies are equal to half the total stored energy N s p of 
the basis function. The volume integral involving 
is transformed to a surface integral thanks to the iden- 



tity V 



(V^) • E* g + 



,V • E; //S . The 



electric field basis functions are divergence-free, so only 
the surface integral j gv <^E* /(3 • e_|_ dS remains. The 
surface integral over the two cross-sections of the wave- 
guide vanishes because both <j)p and E s /y3 satisfy the Flo- 
quet condition. The lateral metallic parts are at imposed 
time-invariant potentials, a property <f)p inherits, so the time 
derivative of the corresponding surface integral vanishes. 
This concludes our demonstration that the evolution of the 
V s p and I s p representing the propagating field is decoupled 
from the space-charge field 
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We perform the inverse /3-transform of (fTTt and (fT2l . Like 
the usual Fourier coefficients, the /3-transform of a product 
is the convolution of its factor transforms. The source term 
in Maxwell- Ampere is transformed by © into an integral 
over the complete volume V% : 

/ J,-E: /3 d 3 r= f J2 J (* + nde *)- E U*)j nl3d d 3 r 

= / J(r + nde x ) ■ B* sf! (r + nde x ) d 3 r 
= / J-E:,d 3 r, 
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where the second equality follows from the Floquet prop- 
erty of E s/ 3. Finally, using the inverse /3-transform, 
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J-F s ,_„d J r, (13) 
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where 



' r E -^e inf3d d(pd) (15) 
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whose /3-transform is F S(3 = E^/iV,^. 
The electric and magnetic fields, given by 

E(r,t) = Vy, B (f)E,,_ n (r)- V^(r,t), (16) 



s,n 



H(r,t) =i ^J m (t)H s ,_ n (r), 



(17) 



are rewritten 

E(r,t) = £V m (t)E s0 (r - ncfe x ) - V0(r, t),(18) 



H(r,t) = i^I s „(t)H s0 (r - nde x ). 



(19) 



Each eigenfunction H s/3 derives from a vector potential 



/i H s/3 = rot A s/9 . 



(20) 



The inverse /3-transform yields /i H sn = rot A sn . Using 
(0, we find 



A = i^/ s „A S! _„ + Vu 



(21) 



where u is any function of space corresponding to a gauge 
choice. 

These expressions offer an interpretation for the index 
n. The electric field at a particular position is the super- 
position of the field shape E s0 generated by the adjacent 
cells of the periodic structure with amplitudes V sn . There- 
fore V sn (resp. I sn ) can be seen as the contribution of the 
s-th eigenmode in cell n to the electric (resp. magnetic) 
field. The propagation is described by linearly coupled har- 
monic oscillators. Cell n is coupled with cell m through 
the coefficient ft s „ = £ U.p e i( - n - m ^ d d{fid), 
viz. the corresponding Fourier coefficient of the disper- 
sion curve O s/ 3. The reciprocity condition Q Si -p = Q s/ 3 
ensures f2 Si _„ = O s „ E R + . The same condition 
on the fields reads E s ._^ = E*^, therefore E s „ = 

±J+:W s ,j^d(pd) = I» (/^ E^e^-dCScO 
is a vector with real components only. Symmetrically, 
iH sn has only real components. In summary, V sn , I sn , 
E sn , F sn , iH s „, iA s „, U sn , N sn are all real quantities. 



Starting from this model, we can now construct a hamil- 
tonian for the electromagnetic field. The standard 
hamiltonian \ ^(eolEp + /z |H| 2 )d 3 r results from the 
Lorentz force and the Maxwell-Ampere equation. In this 
work, we want to make no assumption on the final form of 
the electromagnetic hamiltonian because corrective terms 
might be needed. So we cannot use it as a starting point. 
Instead, the previous development provides a hamiltonian 
for the fields alone, independently from how they are cou- 
pled with the electrons. The vectors V s = (. . . V sn . . . ) 
and I s = (. . . I sn . . . ) represent the state of the fields, 
and their time evolutions derive from the hamiltonian (with 
conjugate variables V, I) 

n em (V 8 , QseN = \ J2(V 8 Q S V* + I s Q s lt), (22) 



where •* denotes vector transpose, and Q s is the infinite 
matrix with entries {Q s ) n . m = Q s . n - m for (n,m) £ Z 2 . 
Till now, the basis functions E s/3 and H s/ 3 can have arbi- 
trary physical dimensions. But adopting (l22l as a hamilto- 
nian results in V S Q S V* and I s Q s Il being energies. Hence 
V s and I s have the same dimensions. Equations (fT3l-(fT4l 
impose £l sn to be a frequency. Therefore V s and / s have 
the dimension of the square root of an action (energy di- 
vided by frequency). The electric and magnetic field have 
their standard dimensions and N s p has the dimension of a 
frequency. 

The hamiltonian of the electron in the field of a given 
4-potential (0, A) is ifTlll 



%ci(p,r) 



m 2 c 4 + c 2 |p 



eA(r)| 2 + e(f)(r) (23) 



where e and m are the electron charge and mass, c is the 
celerity of light, and the momentum p is the canonical con- 
jugate of the electron position r. The complete hamiltonian 
for the fields and M electrons is necessarily the sum of the 
purely electromagnetic hamiltonian H cm with the Lorentz 
force dynamic hamiltonian Tl e \, where A is a function of 
the I sn as given by d2TI) and generates the space-charge 
field, 

^(Pjj v ji V s , Is){l<j<M,seN} = %el( r j> Pj") + %cm- 

(24) 

The complete physics of the system is governed by it, in- 
cluding the source term in the Maxwell-Ampere equation. 
This source term has not been used to build the hamiltonian 
and is equivalent to equation ( fl"3l) . In presence of a single 
electron with position r(i) and velocity v(i), we have fill 
§29] J(r, t) = ev(t) 5(r - r(t)), and (O simplifies into 



V„ 



-ev • F s _„(r) - ^2 VL sm L 



S.n — 771' 



(25) 



Its derivation from % is obtained considering that V sn is the 
generalized momentum canonically conjugate to the gener- 
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alized coordinate I sn , 
dH 



dlsn 



+ 



eA 



dl sn ' ^m 2 c A + c 2 |p - eA| 2 



■ eiA s 



(26) 



: r = dH cl /dp 

eA| 2 , that derives directly from 



Here we used the relation v 
eAJ/^mV + c 2 |p 
the electron hamiltonian (l23l l. 

Equations d25l) and d26l) are incompatible unless 

F sn (r) = -iA sn (r). 

The /3-transform gives E s ^/iV s ^ 



c 2 (p 



definition, E s/3 = - 
tions reconcile when 



(27) 

-iA s/3 . By 
Vu S( g, so the two equa- 



s/3 



N. 
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ft 



and Vu 



s/3 



0. 



(28) 



It is always possible to choose basis functions satis- 
fying the Coulomb gauge, which is the second condi- 
tion. The first condition clarifies the question of the 
electromagnetic energy. In this work, it is given by 
the hamiltonian H em while the usual Poynting energy 
is | L (e |E| 2 + /io|H| 2 )d 3 r. In the latter expres- 
sion, consider the electric part £ = |e j v |E| 2 d 3 r = 
^oY^nlvo l"^( r n de x )\ 2 d 3 r. Using Parseval's rela- 
tion, this is also |e j v (E^ 2 d(/3d) d 3 r. Now, us- 
ing the decomposition of the electric field, we find 

£ = l e ° Jvo C I E s Vs^sp? d(Pd) d 3 r. Orthogonal- 
ity of the eigenmodes reduces this expression to £ = 
\ S s fa* N s p\V s /3 1 2 d(/3d). A similar calculation applies 
to the magnetic energy. Finally, using the inverse f3- 



transform and the condition N. 



s/3 



f2 S(3 , we find that both 
expressions for the energy, "H cm and Poynting's, are equal : 

\ J (e |E| 2 +/, |H| 2 ) d 3 r = \ Y]y*Q.Vl+LQA 

(29) 

In this work, we constructed a hamiltonian describing 
the electrons coupled with the electromagnetic fields prop- 
agating in a periodic structure. This coupled oscillator 
model may be viewed as a generalization of the telegraph 
delay line where each period consists of one inductor and 
capacitor, coupled with their nearest neighbours and elec- 
trons lfl2ll . In contrast, in the present work all periods are 
coupled with each other. This model enables computing the 
time-dependent behaviour of a helix traveling-wave tube, 
faster than industrial PIC codes 

While we focused on divergence-free, propagating 
fields, the space charge field may be taken into account by 
adding the "space charge hamiltonian" of two or more in- 
teracting electrons or by adding Darwin's first-order rel- 
ativistic approximation [8, §12.6] (Green functions are 



needed to take into account the waveguide boundary condi- 
tion). In this action-at-a-distance approach with two elec- 
trons interacting through their electrostatic field, the prob- 
lem of infinite electromagnetic mass does not occur. This 
hamiltonian is formally independent of time but depends 
explicitly on space through the A sn (r). Indeed the bound- 
ary conditions - for example a corrugated metallic wall - 
enable the field to exchange momentum but not energy. 
Maxwell's equations and the Lorentz force correctly ac- 
count for the recoil from the radiated field. This hamilto- 
nian is the sum of two unconditionally positive contribu- 
tions which are minimum when the system is at rest, a situ- 
ation precluding self-acceleration. As a consequence of the 
hamiltonian, the source term in the equation of Maxwell- 
Ampere and the Lorentz force consistently express how the 
electron is coupled with the propagating field. 
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